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Chalcogenide alloys are materials of interest for optical recording and non-volatile memories. We 
perforin ab-initio molecular dynamics simulations aiming at shading light onto the structure of 
amorphous Ge2Sb2Tes (GST), the prototypical material in this class. First principles simulations 
show that amorphous GST obtained by quenching from the liquid phase displays two types of short 
range order. One third of Ge atoms are in a tetrahedral environment while the remaining Ge, Sb 
and Te atoms display a defective octahedral environment, reminiscent of cubic crystalline GST. 



Phase change materials based on chalcogenide al- 
loys are presently used in optical storage devices 
(DVD) and are promising materials for non-volatile 
electronic memories^. Both applications rely on the 
reversible and fast transition between the amorphous 
and crystalline phases which have different optical 
and electronic properties. Among the chalcogenide 
glasses, Ge2Sb2Tc5 (GST) is the material of choice 
for non-volatile memory applications due to its supe- 
rior performances in terms of speed of transformation 
and stability of the amorphous phase 1 . 

In spite of its great technological importance, the 
microscopic structures of amorphous GST (a-GST) 
and the detailed mechanism of the phase transforma- 
tion are largely unknown 2 ". In the past the structure 
of a-GST was implicitly assumed to be a disordered 
version of the metastable cubic (rocksalt) crystalline 
geometry. Very recently based on EXAFS/XANES 
measurements, Kolobov et al£ proposed that in a- 
GST Ge is tetrahedrally coordinated as opposed to 
its octahedral coordination in the crystalline phases 
(hexagonal and metastable cubic) . Based on ab-initio 
calculations Welnic et aJ^ proposed a spinel-like ge- 
ometry for the local structure of a-GST. However, the 
model of a-GST proposed by Kolobov et al is in con- 
trast with other interpretation of EXAFS data^ and 
with more recent Reverse Monte-Carlo (RMC) models 
fitted to XRD data£. On the other hand, models pro- 
duced by RMC are subject to large uncertainties when 
fitted to the total scattering function only, which for 
a-GST is the weighted sum of six partial pair correla- 
tion functions. Further investigations are required to 
obtain a more compelling characterization of a-GST 
which would facilitate the search for better performing 
materials. 

In this respect, first principles atomistic simulations 
can provide precious insight. On the basis of success- 
ful previous studies on the crystalline phases^^ it is 
to be expected that Density Functional Theory works 
well also for the amorphous GST investigated here. 

Ab-initio Molecular Dynamics simulations have 
been performed using the scheme of Kiihne et a^. 
In the spirit of the Car-Parrinello (CP) approach the 
wavefunction is not self-consistently optimized during 
the dynamics. However, in contrast to CP, large inte- 
gration time steps can be used in the simulation. This 



scheme leads to a slightly dissipative dynamics of the 
type — 7dR/ where R/ are the ionic coordinates. In 
Ref. 10 it is shown how to compensate for this dissi- 
pation and obtain a correct canonical sampling. This 
scheme has been implemented in the CP2K suite of 
program s 11 ! 12 . We use the PBE exchange correlation 
functional^ and Goedecker-type pseudopotentials^. 
The Kohn-Sham orbitals are expanded in a TZVP 
Gaussian-type basis set and the charge density is ex- 
panded in a planewave basis set with a cut-off of 
100 Ry. 

The initial configuration is the metastable cubic 
GST where Te occupies one sublattice of the rock- 
salt crystal and Ge, Sb and vacancies are randomly 
placed in the other sublattice in an orthorhombic su- 
percell of size 21.97 x 21.97 x 18.63 A 3 (270 atoms) 
at the density of 0.030 atoms /A 3 close to the ex- 
perimental value for a-GST—. The system has then 
been heated and equilibrated for 6 ps at 2300 K and 
then quenched in 16 ps and further equilibrated for 18 
ps at 990 K. The parameter 7£> = 4 ■ 10~ 4 fs^ 1 has 
been determined as in the Ref. OU The static proper- 
ties of liquid GST turned out to be in good agreement 
with experiments^^ as detailed in the supplementary 
materials (Figs. S1-S3 in EPAPS (Ref. M))- In or- 
der to generate a model of a-GST the liquid has been 
brought to 300 K in 18 ps. The calculated x-ray scat- 
tering function S(Q) of a-GST reported in Fig. Q] is in 
good agreement with XRD data£, especially consider- 
ing that the latter are on the as-deposited amorphous 
film which might be structurally slightly different from 
a-GST quenched for the melt as suggested by their 
different optical reflectivity^. Incidentally, we recall 
that a-GST of relevance for applications is the phase 
quenched from the melt for which, to our knowledge, 
no XRD data are available. Average coordination 
numbers for the different species of a-GST (Table I) 
are computed from partial pair correlation functions 
(Fig. S4 in EPAPS (Ref.©). The distribution of the 
coordination number for the different species are re- 
ported in the inset of Fig. [2] Ge and Sb atoms are 
mostly 4-coordinated and form bonds preferentially 
with Te atoms. However, we observe a large fraction 
of homopolar Ge-Ge, Sb-Sb and Ge-Sb bonds, namely 
38 % of Ge are bonded with at least another Ge or 
Sb, a percentage which raises to 52 % for Ge atoms 
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FIG. 1: (Color online) Calculated and experimental x-ray 
scattering factor S(Q) of a-GST. The theoretical S(Q) has 
been computed from the partial structure factors Sij(Q) 
weighted by the Q— dependent x-ray atomic form factors^ 
(S(Q) - 1 = F X (Q) of Eq. 56 in Ref. \M)- Sij(Q) are 
obtained in turn by Fourier transforming the partial pair 
correlation functions (3 ps at 300 K). Inset: a sketch of the 
geometry of defective octahedral sites of 4-coordinated Ge 
and 3-coordinated Te. 



4-coordinated. The fraction of Sb bonded to another 
Sb or Ge is instead 59 %. Homopolar Ge-Ge bonds 
have been recently detected by x-ray fluorescence in 
cubic GST grown epitaxially on crystalline GaSb^. 
The concentration of Te-Te bonds is somehow lower, 
27% of Te atoms are involved in homopolar Te-Te 
bonds arranged into dimcrs and trimers (see Fig. S5 
in EPAPS (Ref. [H)). The large concentration of ho- 
mopolar bonds is not reproduced by the RMC model 
— . Insight on the local geometry is further gained 
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total 


Ge 


0.275 


0.270 


3.277 


3.823 


Sb 


0.270 


0.588 


3.166 


4.025 


Te 


1.311 


1.267 


0.288 


2.866 



TABLE I: Average coordination number for different pairs 
of atoms computed from the partial pair correlation func- 
tions (Fig. S4 in EPAPS (Ref. 01)). 

from the angle distribution function in Fig. [21 The 
broad peak at ~90° and the weaker structure around 
~170° are reminiscent of the distorted octahedral- like 
geometry of the metastable cubic crystal. A snapshot 
of the a-GST model is shown in Fig. [3] Angles at 
^90° and ^180° clearly dominate the bonding net- 
work. For Te, only angles at ^90° are found (see 
inset of Fig. [I}. However, the main coordination of 



Ge(Sb) and Te of four and three, respectively, is lowe? 
than the ideal octahedral value of six. In our model 
the presence of neighboring vacancies is responsible 
for the lower coordination while the bonding angles 
remain close to ^90° and ~180° as in the metastable 
cubic phase. Te and Sb atoms adopt this configura- 
tion, while only a fraction of Ge atoms are in a defec- 
tive octahedral site. A large fraction of 4-coordinated 
Ge atoms are in a tetrahedral environment as inferred 
from EXAFS/XANES measurements^. A signature of 
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FIG. 2: (Color online) Angle distribution function (to- 
tal and resolved for different central atoms). Inset: dis- 
tribution of coordination numbers of different species ob- 
tained by integration of the partial pair correlation func- 
tions (Fig. S4 in EPAPS (Ref.GJ)). 

the tetrahedral geometry is already visible as a shoul- 
der at ^110° in the angle distribution function for Ge 
atoms in Fig. [H A better indicator of the tetrahe- 
dral geometry is given by the local order parameter— 
q = 1 — |X)i>fe(i + cosOijk) 2 where the sum runs 
over the couples of atoms bonded to a central atom j. 
The distribution of the local order parameter q for Ge 
atoms is reported in Fig. [3] for different coordination 
numbers. The 4-coordinated Ge distribution is clearly 
bimodal with peaks corresponding to defective octa- 
hedra and tetrahedra. In contrast, the ^-distribution 
for Te and Sb does not show any signature of the tetra- 
hedral geometry (Fig. S6 in EPAPS (Ref. HI)). The 
^-distribution for 4-coordinated Ge is further analyzed 
in terms of atoms bonded to Te only or to at least one 
Ge or Sb. The presence of bonds with Ge or Sb clearly 
favors the tetrahedral geometry. Only few tetrahedral 
Ge are bonded to Te only. On the other hand, all 
Ge with more than one homopolar bond (with Ge or 
Sb) are in the tetrahedral geometry. By integrating 
the tetrahedral peak of the g-distribution in the range 
0.8-1.0 we estimate that 33 % of Ge atoms are in a 
tetrahedral environment. The average bond length of 
Ge is slightly shorter in tetrahedral sites (2.73 A) than 
in defective octahedral sites (2.82 A, see Fig. S6 in 
EPAPS (Ref. 16)). The presence of tetrahedral sites 
and the absence of 6-coordinated octahedral sites (or 
distorted octahedra sites with 3+3 coordination as in 
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FIG. 3: (Color online) Distribution of the local order pa- 
rameter q for Ge (see text). q=l for the ideal tetrahe- 
dral geometry, q—0 for the 6-coordinated octahedral site, 
and g=5/8 for a 4-coordinated defective octahedral sites 
Top panel: g-distribution resolved for Ge with different 
coordination number. Bottom panel: g-distribution for 4- 
coordinated Ge further resolved for Ge with at least one 
homopolar bond (with Ge or Sb) or bonding with Te only 
(no homopolar bonds). Top panel inset: ring distribution 
function of a-GST computed as in Ref. 20i. Bottom panel 
inset: snapshot of the 270-atoms model of a-GST. 



crystalline a-GeTe2i) is consistent with the model in- 
ferred from EXAFS/XANES data in Ref. 3. The co- 
existence of tetrahedral and defective octahedral sites 
found here might reconcile the different interpreta- 
tion of EXAFS data^ based on the assumption of 



a unique environment for all Ge atoms. The inclusion 1 
of both configurations in a revised fitting might pro- 
vide a better agreement with experiments. Turning 
now to the medium range order of a-GST, we report 
the ring distribution in Fig. [3J Both even and odd 
rings are present, as opposed to the RMC results — . 
The distribution has a pronounced maximum at four- 
membered rings typical of the rocksalt structure. The 
presence of very large rings (with very large aspect ra- 
tio) reveals a somehow more open and less connected 
structure than that of the cubic crystalline phase. 

In summary, based on first principles simulations we 
have provided novel insight into the structure of amor- 
phous Ge2Sb2Te5 quenched from the melt. Most of Ge 
and Sb atoms in a-GST are 4-coordinated while Te is 
mostly 3-coordinated in defective octahedral-like sites 
which recall the local environment of cubic crystalline 
GST. However, as many as 33 % of Ge atoms are in a 
tetrahedral geometry, absent in the crystalline phase, 
and favored by the presence of homopolar (Ge-Ge and 
Ge-Sb) bonds. The coexistence of the two types of lo- 
cal environment is the key to understand the two ap- 
parently contradictory and peculiar features of GST 
exploited in the devices, namely the strong optical 
(and perhaps electronic) contrast between the amor- 
phous and crystalline structures and the high speed of 
the phase change. This conclusion is also supported 
by the recent theoretical work in Ref. 0) on models of 
amorphous GeSbTe alloys showing the presence of a 
fraction of Ge in tetrahedral sites is indeed sufficient to 
produce a strong optical contrast between the amor- 
phous and crystallie phases. 

During the preparation of the manuscript we have 
received a preprint by Akola and Jones where simula- 
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their results are similar to ours. Computational re- 
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Parallelo 2007". S.C. acknowledges financial support 
given by Fondo Sociale Europeo, Ministero del Lavoro 
e della Previdenza Sociale and Regione Lombardia. 
Discussion with R. Bez, A. Modelli, A. Pirovano and 
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FIG. SI: Calculated and experimental scattering functions 
S(Q) of liquid GST from x-ray (953 K, Ref. 6) and neutron 
diffraction (1073 K, Ref. 15). S(Q) has been calculated ac- 
cording to the definition by Keen (Ref. 23). Namely, S(Q) 
is the sum of partial structure factors Sij(Q) weighted by 
Q— dependent atomic form factors (S(Q) — 1 = F x (Q) of 
Eq. 56 in Ref. 23) or by neutron scattering lengths (Eq. 19 
in Ref. 23). Sij(Q) are determined in turn by Fourier 
transforming the partial pair correlation functions gij(r) 
(see Fig. S2, 990 K). 
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FIG. S2: Total and partial pair correlation functions of liq- 
uid GST (18 ps at 990 K). The first peak in Sb-Sb and Ge- 
Ge partial pair correlation functions are due to homopolar 
bonds. We have checked that the system was liquid from 
the linear asymptotic growth of the mean square displace- 
ment. From the Einstein relation we found D — 4.88- 10 
cm 2 /s. This is close to the estimate of D — 4.55 • 10 -5 
cm 2 /s obtained from shorter standard Born-Oppenheimer 
calculations performed to check the validity of the 
choice. To our knowledge no experimental data on the 
self-diffusion coefficient of liquid GST are yet available. 
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FIG. S3: Angle distribution functions of liquid GST at 990 
K (total and resolved for different central atoms). Angles 
at 60° (absent in crystalline and amorphous GST) corre- 
spond to triangular rings. 
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FIG. S4: Total and partial pair correlation functions of 
a-GST (3 ps at 300 K). The vertical lines are the bond- 
ing cutoff used to define the coordination numbers (Table 
I and inset of Fig. 2). For Sb-Te the bonding cutoff has 
been chosen as the outer edge of the partial pair correlation 
function of our model of cubic GST (300 K). The maxi- 
mum of the GeTe and SbTe pair correlation functions are 
2.79 A and 2.94 A to be compared with the experimental 
values of 2.61 A and 2.85 A, respectively (EXAFS, Ref. 3). 
The first peak in Sb-Sb and Ge-Ge partial pair correlation 
functions is due to homopolar bonds. We have checked 
that the cutoff distances for Sb-Sb and Ge-Ge correspond 
indeed to the formation of covalent bonds by looking at the 
position of the Wannier orbitals (N. Marzari and D. Van- 
derbilt, Phys. Rev. B 56, 12847 (1997)). The second 
peak in Sb-Sb and Te-Te correlation functions correspond 
to the diagonal of the ABAB square ring of the rocksalt 
geometry. The outer peaks above 6 A correspond to the 
diagonal of two planar, edge sharing ABAB squares. In- 
stead, the peak around 5 A expected for the diagonal of 
the elemental rocksalt cube is lacking in Ge-Te and Sb-Te 
correlation functions. 
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FIG. S7: Total pair correlation function for Ge in tetrahe- 
dral sites (continuous line) or in defective octahedral sites 
(dashed line). Ge-X bonds are shorter for Ge in tetrahe- 
dral environment. 



FIG. S5: Snapshot of a-GST showing Te (black spheres) 
and Sb (gray spheres) atoms in dimers and short chains. 
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FIG. S6: Distribution of the local order parameter q (see 
text) for Sb and Te in a-GST. Vertical lines indicate the 
values of q for selected ideal geometries. No Te and Sb 
atoms are in tetrahedral environments. 



